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Abstract 

We examine statistical field theories of polymeric fluids in view of performing numerical simula- 
tions. The partition function of these systems can be expressed as a functional integral over real 
density fields. The introduction of density field variables serves to decouple interactions among 
non-bonded monomers, and renders the resulting effective Hamiltonian H for the field theory real 
and the Boltzmann factor exp(— H) positive definite. This is in contrast with conventional (Ed- 
wards type) field theories expressed in terms of chemical potentials that have complex H. The 
density field theory involves the calculation of an intermediate functional integral, which provides 
the entropy of the polymer fluid at a fixed density profile. This functional integral is positive 
definite and we show that in the thermodynamic limit of large systems, it is dominated by saddle 
points of the integrand. This procedure side-steps the "sign problem" in the chemical potential 
field formulation. The formalism is illustrated in the context of models of flexible polymers. We 
discuss the implications for field-theoretic computer simulations of polymeric fluids. 
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I. INTRODUCTION 



Field theories have proved to be a useful theoretical framework for studying the equi- 
librium and time-dependent properties of a wide variety of complex fluid systems. In the 
1960s, Edwards introduced field-theoretic methods to the field of polymer physics and ap- 
plied these techniques fruitfully to examine the universal aspects of chain conformational 
statistics and thermodynamic properties of polymer solutions Q]. Perturbation methods, 
mean-field approximations, and variational schemes were introduced for tackling important 
problems such as the excluded volume effect and addressing collective phenomena, such as 
screening in semi-dilute solutions. These analytical methods were considerably enhanced in 
the 1970s and 1980s with the advent of renormalization group theory, the introduction of 
scaling concepts, and the discovery of the connection of polymer statistical mechanics with 
the n- vector field theory model for n — > Q]. 

Unfortunately, most of these analytical techniques are difficult to apply in the important 
case of inhomogeneous polymers, which include multicomponent polymer alloys, block and 
graft copolymers, nano-composites, thin films, emulsions, and suspensions. Appropriate 
field theory models for such systems are straightforward to construct, but the analytical 
tools available for examining the properties of the models are rather limited. One powerful 
tool is self-consistent field theory (SCFT), which amounts to a saddle point approximation 
to the partition function of the field theory 0, 0]. However, the SCFT field equations 
determining the "mean fields" are nonlinear, nonlocal, and defy analytical solution. Certain 
limiting cases of SCFT, such as weak inhomogeneities or fields that vary slowly over the 
size of a polymer, and for simple geometries, are analytically tractable. In such special 
circumstances, however, the SCFT solutions still reflect the mean-field assumption intrinsic 
to the theory, which further restricts their applicability. In particular, such solutions are 
primarily limited to the description of concentrated solutions or melts of high molecular 
weight polymers. Even in the case of dense molten blends or copolymers, SCFT can fail in 
situations where soft modes appear, such as near critical points, unbinding transitions, or 
in micellar or microemulsion phases Q], since thermally- induced field fluctuations can then 
be significant. 

In recent years, it has become apparent that field theory models of polymers and complex 
fluids, besides serving as the platform for analytical calculations, can also provide a flexible 
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framework for numerical simulations |6j. Building on significant advances in numerical 
methods for solving the mean-field SCFT equations R, 0], techniques have been developed 
for direct numerical sampling of field fluctuations |(| |J. These methods open many exciting 
possibilities for carrying out "field-theoretic simulations" of a wide variety of inhomogeneous 
polymer systems. Aside from the numerical errors associated with representing the fields 
and the usual statistical errors intrinsic to any computer simulation technique, field-theoretic 
simulations provide a route to studying the exact properties of a polymer fluid model in 
the absence of any simplifying approximations. Such field-based simulations also facilitate 
systematic coarse-graining, making it possible to equilibrate much larger systems with a 
broader range of structural scales than in conventional atomistic computer simulations. At 
the present time the methods are limited to equilibrium structure and properties, although 
research is underway to extend them to nonequilibrium situations. 

A difficulty with the current strategies for field-based simulations of polymeric fluids 
is that they rely on the introduction of Hubbard-Stratonovich "chemical potential" fields, 
first employed for polymers by Edwards to decouple the non-bonded interactions among 
monomers. This allows one to explicitly trace out the monomer coordinates in favor of 
the introduced fields, generating a field theory model for a polymer solution or melt. Un- 
fortunately, this description leads to an effective Hamiltonian Q{4>), expressing the energy 
of a chemical potential field configuration 0(r), that is complex, rather than real. As a 
consequence, the Boltzmann factor exp(— Q) is not positive definite, so conventional Monte 
Carlo sampling of the partition function is problematic ^^|. This difficulty in sampling field 
configurations is a special case of a more general "sign problem" , which involves the devel- 
opment of techniques for efficient numerical evaluation of high dimensional integrals with 
non-positive definite integrands. While promising new methods have appeared for battling 



the sign problem 
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I12I ] , the problem is by no means "solved" . Oscillations in the phase 
factor exp(— i Im Q) severely limit the convergence of field-based simulations of polymeric 
fluids. 

In the present paper, we describe an alternative formulation of polymer field theory that 
essentially side-steps the sign problem. Rather than retaining just a Hubbard-Stratonovich 
chemical potential field for each monomer species, we also retain a second density field. 
This two field approach is not new J^l Q] , but rather is a usual starting point for deriving 
SCFT. However, we show in the present paper that at a fixed value of the density field, the 
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functional integral over the chemical potential field can be performed in the saddle point 
approximation, and that this approximation is exact in the thermodynamic limit of large 
systems. This we believe to be a new result and one that has profound consequences for 
field-theoretic computer simulations of polymeric fluids. 



II. BASIC FORMALISM: POLYMER SOLUTIONS 
A. Transformation to a Field Theory 

Consider a solution of M homopolymer chains in a volume V, each chain composed of 
N monomers of size a. The solvent is not considered explicitly, the temperature is denoted 
by T, and (3 = 1/T. The formalism presented below can be easily generalized to blends of 
different chains and to melts or solutions of copolymers as will be seen later. 

We assume that the monomers interact through a 2-body interaction potential (of mean 
force) denoted by v(r). In the celebrated Edwards model, the interaction v(r) is modelled 
as a Dirac delta function, S(r), but in general, the interaction v can represent both sol- 
vent effects, electrostatic effects or any other effective monomer-monomer interaction. The 
canonical partition function of the system reads 



(1) 

where f Vr k denotes a path integral over all conformations of the fcth polymer chain. 

Next, we introduce a monomer (number) density field p(r) as a new integration field. 
The microscopic expression for the monomer density, p(r), is given by 

m N 

p{r) = Y J / ds5(r-r k (s)) (2) 
k=i Jo 

If we constrain p(r) to p(r) at each point r of space by the use of a 5-function, the partition 
function of the polymer system can be written as 

Z — J Dp exp^— drdr' p{r)v{r — r')p{r')^ Z{p) (3) 



4 



where the function Z(p) is defined by 

= m I n ^ «p (-^ t jf * (£) ! ) n * (pf) - w) <4) 

and J £>p denotes a functional integral over the real, scalar density field p(r). The object 
Z(p) is the partition function of a set of M continuous Gaussian chains of N monomers, 
constrained in such a way that the total monomer density is fixed to p(r) at any point r 
of space. This partition function evidently counts the total number of chain conformations 
that are consistent with the density profile p(r). We thus write 

Z(p) = exp[S(p)} (5) 

and identify S(p) with the entropy of a system of M chains constrained to have a density 
profile p{r). From its very definition, it is clear that Z(p) is real and positive for any real 
p(r). In addition, we expect Z(p) to scale with V, N, and M as (Vg N ) M /Ml, where V is 
the volume of the system and g is some positive, real constant. 

At this stage, it is clear that we have expressed the partition function of the polymeric 
system in terms of a functional integral over a real field with a positive definite Boltzmann 
weight. We have yet to show how to compute the entropic part of the weight, namely the 
factor Z(p). 

Using a Fourier representation of the delta function, the above expression for Z(p) can 
be written as 



Z(p) = jfi f £W>exp (i J dr^(r)p(r) + M\nQ(i(p) 



(6) 



where 

Q{i4) = I Pre-i^MS)^/****)) (7 ) 



is the partition function of a single polymer in the presence of the potential i<p(r). It is 
important to note that this is a purely imaginary potential because the functional integral 
in eq. (JBJ) is over a real field 4>(r). We recognize in eq. (J7J) a Feynman path integral for 
a quantum mechanical problem with quantum Hamiltonian Ti = -^V 2 + i<fi(r). Using 
standard quantum mechanical notation, we may write 

Q(i(j)) = [ dr [ air' < r\e' NH \r' > (8) 
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B. Expansion about the Saddle Point 



Of particular interest is an asymptotic expansion of Z(p) for large M to examine the 
thermodynamic limit of M — > oo, V — > oo, M/V finite. For this purpose, we perform 
a saddle point expansion of the integral. The first step is to write an equation for the 
mean-field (saddle point) </> ( r ) that produces the desired p(r) 

P(r) = -M S -Mlf < 9 > 
M f N 

' ds V(r,s)V{r,N - s) (10) 



Q(i<Po) Jo 

Here, s is denotes a contour location along the chain and the "propagator" $(r, s) satisfies 
the diffusion equation (analogous to the Schrodinger equation) 

(-Q- s -jV 2 + iMr))V(r,s) = (11) 

subject to the initial condition \l/(r, 0) = 1. The single chain partition function is related to 
the field * by 

Q(i0) = f dr V(r,N) (12) 



^From the structure of the above equations, we can make a few general statements. First, 
a necessary condition for eq. to have a solution 4> (r) is that p(r) should be positive at 
any point r and that its integral over space should be equal to MN. It is also clear from 
eqs. (0) and (fTUj) that a self-consistent solution may exist with a real i<p Q (r), i.e. purely 
imaginary 4> {r). From the structure of these equations, we see that if </> ( r ) ^ s a mean- field 
(saddle point) solution, one can add to it any constant and obtain a new solution. To lift this 
degeneracy of the saddle-point, we may specify the value of the field <ft at a specific point 
in space, or its spatially-averaged value. Having made such a choice, it can be proven that 
such a real solution i(f> Q (r) not only exists, but is a unique solution of the above equations. 

To perform the asymptotic expansion about the saddle point <j) at higher orders, we shift 
the integration variable by writing 

0(r) = o (r) + x(r) 

The partition function then reads 

Z{p) = Z {p) J V X exp (i J dr x(r)p(r) + Mln Q{i X )\ (13) 



where 

[Vre-^Z H^) 2 -* fo N ds Mr(s))~i j" ds X (r(s)) 

and Z denotes the mean-field value of the partition function 

MP) = ^yexp (i j dr<f> {r)p(r) + MhxQ(i<f> )) (15) 

Next, we expand the exponent in the integrand of eq. (JT3jl in powers of \. For this 
purpose, it is convenient to reexpress eq. (fTljl as 

Q(i X ) =< e -ifo N dsx(r(s)) >q ( lg ) 

where the notation < ... > describes the expectation value over the conformations of a 
single polymer subjected to the mean-field potential i<fr . Expanding in powers of x an d 
rescaling the field x by a factor of vM, we obtain 



Z(p) = Z (p) J V X exp (~ J drdr' X (r)G 2 (r,r') X (r') + W( X )) (17) 

where 

W{ X) = E p!MP/2-i J dr ± ■ ■ ■ dr P *( r i) ■ ■ ■ x(r P )G p (r u ...,r p ) (18) 

and G p (r%, . . . , r p ) denotes the p-point connected (cumulant) correlation function for the 
density of a single chain, defined by 

G p (n, . . . , r p ) =< p{r x ) . . . p{r p ) > c (19) 



where p{r) = L ds 5(r — r(s)) is the single chain microscopic density. This expectation 
value is taken with respect to the conformations of a single chain in the mean-field according 
to eq. (jlfij) with the usual definition of connected averages Q]. Note that there is no linear 
term in x appearing in the exponent of eq. ()17j) since the action is stationary with respect 

tO O . 

The partition function of eq. (JT7j) can be calculated in a perturbation expansion in powers 
of 1/M. This is most conveniently done by the use of Feynman diagrams in which propagator 
lines are associated with factors of ^(r, r') = Crg (r, r') and pth-order vertices are assigned 
factors of [{-i) p /M p ^- l ]G p {r u ...,r p ). The result is that 

Z(p) = Z (p) exp ^— -Tr hiG^r, r') + Sum of all Connected Diagrams^ (20) 



where the connected diagrams have one or more vertices of degree p > 3. Since the vertices 
in eq. (j!9|) are connected correlation functions, all terms in the exponent of eq. (j2*Uj) are 
extensive, that is they are proportional to the volume V of the system. The general expansion 
of Z(p) thus takes the form 



Z{p) = exp 



(21) 



where the fj are contributions to the free energy density that are independent of M and V in 
the thermodynamic limit. Since all vertices carry powers of 1/M, the sum of all connected 
graphs is negligible in the thermodynamic limit, and eq. ()20j) yields 

Z{p) ~ Z (p)exp (~TrlnG 2 (r,r')\ (22) 

in the limit M — > oo. We note that G2{r,r') is a functional of p, because it represents the 
2-point connected correlation function of a chain subjected to the mean-field potential i(p 
and O is in turn related to p by the saddle point condition eqs. (|9|)-(fTD|). 

^From the form of eq. (|22jh it would appear that there is a correction to the mean- 
field value Zq(p) of the restricted density partition function Z(p) that is associated with 
Gaussian fluctuations in x( r ) about the mean chemical potential field 0o( r )- Interestingly, 
these corrections from Gaussian field fluctuations play no role in the statistical properties 
of the field p for M —>■ oo. This is demonstrated by two independent routes in the next 
section. 



C. Significance of the Fluctuation Term 

To examine the significance of the fluctuation term |Tr lnG^r, r'), we return to the 
original formulation of the field theory prior to the asymptotic expansion in M. From 
eq. (jHJ) it is clear that we can write an effective, real field theory for the p field as 



Z = JVp exp[-H(p)} (23) 
where the effective Hamiltonian H(p) is given by 

H{p) = ^ J drdr' p{r)v{r - r')p{r') - S(p) (24) 
where S(p) = InZ(p) is given in eq. (0)- One way to simulate such a field theory is through 



a Langevin stochastic dynamics 



15| for p. Note that the density field p(r) is constrained 



to be positive at each point r. A simple way to enforce this condition in the functional 
integral ()23|) is to add an external potential which is infinite when p(r) is negative and 
when p(r) is positive. A practical realization of such a potential can be obtained by adding 
an exponential function to the Hamiltonian, and define 



H\p) = H{p) + J dr e- xp{r) 



where A is a constant, chosen large enough to ensure the positivity constraint. 

If we keep in mind that the total number of monomers is fixed, equal to MN, we may 
use the conserved dynamics (model B) and write 

= TV 2 \ P J dr'v{r - r')p{r', t)- < i<p(r, t) > p -X e - Xp(r >j + r){r, t) (25) 
where we have used 

6Sjpl = 6lo g Z(p) 

S P (r,t) 5p(r,t) <l ^ > " 
and T)(r,t) is a Gaussian, white noise with vanishing first moment, (r](r,t)) = 0, and second 
moment satisfying the usual fluctuation-dissipation theorem for a conserved noise, 

< v (r, t)r)(r' 7 1') >= -2TV 2 5(r - r')5(t - t') (26) 

The Laplacian in front of the r.h.s. of (J25j) ensures that the space integral of p(r,t), which 
is equal to the total number of monomers MN, remains constant in time. In the above 
equation, the expectation value < i(f)(r) > p denotes an equilibrium average of the field i(j)(r) 
at constrained p(r) using the statistical weight contained in eq. © for Z(p). 
With the use of eq. (fT7j). it is easily seen that 

< i<f>(r) > P = i4>o( r ) + -t= < X(r) > P 



where 

< > _ J Vx X(r) exp (-| / drdr'x(r)G 2 (r : r')x(r') + W(x)) 
J T>x exp (— \ J drdr'x(r)G2(r,r')x(r') + W(x)) 
and W(x) is given by eq. (fTHj). From this expression it follows immediately that 

< z0(r) > p = i(f> (r) + O(^) (27) 



and thus, in the Langevin equation formalism, it is sufficient to evaluate < i(j){r) > p at the 
mean-field level to capture the leading asymptotic behavior in the thermodynamic limit. The 
evaluation of the Langevin forces on the right hand side of eq. (J25)) therefore requires only 
the computation of <f) (r,t), which follows from the solution of the saddle point equations 
at prescribed p{r,t). The forces entering the Langevin equation are real (since i(j) is real), 
and thus we have a conventional real Langevin dynamics. The key result is that in the 
thermodynamic limit, there is no signature of the fluctuation term (l/2)Tr lnG 2 (r, r') that 
appears in eq. (|22Jl for Z(p). 

The same result can be directly obtained from eq. ()22)l by means of some straightforward 
manipulations. Indeed, from that equation and eq. ([To]) , we have 

5\nZ(p) 



< i(t>(r ) > P (28) 

, ± . . 15 Tr In G 2 (r,r') 

- 2 Mr.) <29) 



The Trln term can be further simplified. Using some simple identities, we get 

Sp(r ) =J drdr ^r G > (r ' r) 

- J drdrdn 5p(ro) 5{iMri) ^2 (r,r) (30) 



^From our previous definitions, we see that 

5G 2 (r, r') _ 
5(i0 o (ri)) 

and the saddle point condition implies 



-G 3 (ry, n ) (31) 



Wo( r i)) _ 



Spi r ) ~M G *^ < 32 » 
Substitution of these results into eq. (|3T)j) leads to a final expression 

S p ( rQ ) = ^( r °) ~ 2M J rfrrfr ' (ir i G '2 1 ( r i ) r o)G'3(r,r',ri)G2 1 (r',r) (33) 

Again, since G 2 and G 3 are connected Green's functions, the triple integral in eq. (|33[) is 
finite in the thermodynamic limit (non extensive) and this shows that to leading order in 
1/M, the expectation value of <ft is given by its mean-field value O , in agreement with 
eq. (|2T|). In addition, it is easy to see that the correction term calculated here is exactly 
the same as the one that would have been obtained as the leading term of the perturbation 
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expansion using eq. (JT7J). We conclude that the fluctuation term (l/2)TrlnG2(r, r') can 
be safely omitted from the effective Hamiltonian H(p) when computing the Langevin force 
5H/5p for the purpose of simulating large systems. Moreover, it follows that the fluctuation 
correction can also be dropped from the Hamiltonian itself, if Monte Carlo simulations are 
to be employed. 

We summarize with the statement that in the thermodynamic limit of large systems, the 
field theory expressed in density variables can be written as 

Z = Jvp exp[-H'(p)] (34) 

with 

H'(p) = ^ J drdr' p{r)v{r - r')p(r') + J dr e" Ap{r) - / dr i<f> (r)p(r) - M\nQ(i<p ) (35) 
and where the real field i<f) (r) is obtained at a particular p(r) by the solution of eq. 



D. Langevin Dynamics and Monte Carlo: Implementation 



The Langevin equation (j25jl can be discretized both in time and space. A variety of time 
integration schemes for such nonlinear stochastic differential equations are available 

It is of interest to compare this real Langevin scheme with the complex Langevin scheme 
that has been used to numerically simulate field theories of polymer solutions and melts 
expressed in chemical potential [i.e. 4>(r)] variables on . Because the effective Hamiltonian 
Q{4>) in such field theories is complex, one must address the sign problem. A useful strategy, 



originally devised by Parisi and Klauder 19, 20], has been to extend the field variables to 
the complex plane (<f) = <p R + itfij) and write a Langevin dynamics 

M^ = - TB *J§J) + "(''•" < 36) 

y^=- Ti ™wk) (3?) 

where rj(r,t) is a real Gaussian, white noise with the following covariance 

< r](r, t)rj(r', f) >= 2T5(r - r')5(t - t') (38) 



This procedure, if it converges [12J, |2l|, |22j, |23( , will properly describe the stationary fluctu 



ation spectrum of a model with complex Q(<f>). We have found it to be reasonably effective 
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at suppressing phase oscillations in simulations of polymer solutions and melts, although 
sign problem is still very pronounced in dilute and semidilute systems. An advantage of 
eqs. ()36| ) -([37 p is that the computation of the complex force 5Q/5(p requires only a single pass 
at solving the diffusion equation (jllj) (with <p —>■ <p), while computation of STi'/Sp in the 
real Langevin scheme of eq. (J25|) requires multiple passes at solving the diffusion equation 
in order to establish the field O consistent with the current p. Whether the extra compu- 
tational burden of the force evaluation in the present density-based real Langevin scheme 
is offset by avoidance of the sign problem remains to be seen and, indeed, may prove to be 
system and problem dependent. 

Note that in a Monte Carlo implementation, the minor complications due to the positivity 
of p(r) and the conservation of the integral J drp(r) are easy to overcome. First, one 
discretizes space. To implement the conservation of the total number of monomers, one starts 
with a configuration of p(r) that has the right number of monomers, and then one chooses 
randomly two lattice sites, say r\ and T2 ■ The elementary Monte Carlo variation of p(r) 
consists in modifying simultaneously p at two points r\ and r2 according to p{r{) — > p(r\)+5p 
and p(r2) — > p{t?) — Sp . This procedure obviously conserves the total integral of p(r) in 
space. Also, to enforce the positivity of p(r), one just has to reject any variation of p that 
produces a negative p. 

III. MORE COMPLICATED SYSTEMS 

The concepts and methods presented above can be extended in a straightforward way 
to virtually any type of inhomogeneous polymer system. We illustrate in the context of a 
solution of AB diblock copolymers. 

Let / be the volume fraction of A monomers and 1 — / the volume fraction of B monomers 
on each diblock copolymer chain with a total of N monomers. Denoting by vaa{t), vbb{t), 
and Vab (r) the respective potentials of mean force (mediated by the solvent) between the 
different species of monomers, and using the same notations as in the previous section, we 
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may write 



Z = W\ J VpAVpBexp / drdr 'PA( r ) v AA(r - r')p A {r') 
x exp J drdr' p B {r)v BB {r - r')p B (r')J 

x exp y-(3 J drdr'p A {r)v AB {r -r')p B {r'yj Z A {p A )Z B {p B ) (39) 



lwhere the functions Z A (p A ) an d Zb(Pb) are defined by 



Za( Pa ) = /n^exp i^ipEl ds (^) 



and 



z B (p b ) = /n^- P (-Ag^ ds g) 

/ M N \ 

IlM^W-E / <k*(r-r fc (a)) (41) 



X 

It is clear from their definitions that both Z A and Z B are positive definite, and using a 
Fourier representation for the delta functionals, one can again show that in the limit of an 
infinite number of chains (with infinite volume and finite concentration), the corresponding 
"chemical potential" field integrals can be asymptotically evaluated by the saddle point 
method. The resulting expression for the entropy functional 

S(p A ,p B ) = \n[Z A (p A )Z B (p B )] (42) 



is 



s (Pa,Pb) = J dr (i(f) A p A + i(f) B p B ) + MlnQ(i(f) A ,i(j) B ) (43) 

where Q(icj) A ,i(j) B ) = f dr ty(r,N). The propagator ty(r, s) satisfies (again, subject to 
*(r,0) = 1) 

(^- s -jV 2 + w(r,s))^(r,s) = (44) 

where w(r,s) = i(f> A (r) for s G (0,Nf) and w(r,s) = i(j) B (r) for s G (Nf,N). The density 
fields are related to the saddle point chemical potentials, 4> A and <p B , by equations analogous 
to eq. dmj) 

M f Nf 
0(%,#b) Jo 
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M ' N 



QWaMb) JNf 

A second propagator \I/*(r, s) enters these expressions due to the head-to-tail asymmetry of 
diblock copolymers. It satisfies 

{ ~d~s ~ T v ' + w * (r ' s))r * (r ' s) = (47) 

subject to \I/*(r, 0) = 1, where w*(r, s) = z0s(r) f° r s e (0, A"(l — /)) and io*(r, s) = ^^(r) 
for s G (iV(l - /),JV). 

A formalism very similar to this was employed recently, without formal justification, for 



a real Langevin dynamics study of spinodal decomposition in binary polymer blends 
and in studying microemulsion formation in ternary copolymer alloys |25l| . The present 
analysis shows that this is indeed a rigorous approach to simulating fluctuating systems of 
macroscopic size. 

IV. SUMMARY AND DISCUSSION 

In the present paper we have examined statistical field theories of polymeric fluids with 
the objective of identifying field-based models that are well-suited to numerical simulation. 
It has been known for some time that essentially any microscopic (particle-based) model 
of a polymeric fluid can be converted to a field theory in density and chemical potential 
variables 0, Q, e.g. for a one-component system 

Z = Jvp Jvcf) exp[-W(p,0)] (48) 

where the effective Hamiltonian 7i(p, 4>) is a complex functional of the two fields. In the case 
of pair interactions among monomers, T~C(p, <fi) is a quadratic form in p, so the density field 
can be integrated out exactly (for potentials with an inverse) to obtain a statistical field 
theory in <p 

Z = Jv(f> exp[-£(0)] (49) 

In this case, the effective Hamiltonian Q{<f) is complex. This type of field theory has been 
studied using approximate analytical methods for many vears jl| and has recently been 
examined with complex Langevin simulation techniques jy, |9j, [l8| . Unfortunately, the non- 
positive definite character of the Boltzmann factor exp(— Q) gives the integrand an oscillatory 
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character ("sign problem"), which can dramatically hinder the convergence of numerical 
simulations. 

Rather than transforming to a field theory in the single field 0, a more typical strategy for 
studying inhomogeneous polymers is to invoke a mean-field approximation. This amounts 
to evaluating both integrals in eq. P8|) at leading order by the saddle point method. This 
reproduces the familiar self-consistent field theory (SCFT) P and amounts to 

Z^exp[-H(p*,4>l] (50) 

where p* and 0* are the "self-consistent" density and chemical potential fields that satisfy 
5H/5p* = 0, 5H/5(p* = 0. This is a powerful approach, but it neglects all field fluctuations, 
so the utility of the SCFT is limited to highly concentrated systems of high molecular weight 
polymers, far from critical points and phase transitions where soft fluctuation modes occur. 

In the present paper, we have discussed a third approach in which the chemical potential 
field is integrated out of eq. (f4*H|) . leading to a field theory solely in the monomer density 

Z = jvp exp[-F(p)] (51) 

This approach has the benefit of producing a real field theory in which exp(— H) is positive 
definite so that the sign problem is avoided and standard Monte Carlo and real Langevin 
simulation methods can be applied. In the past it was assumed that the integral necessary 
to reduce eq. (pIHj) to eq. (J5T|) was intractable, leading various authors to impose a variety 



of approximations in order to construct the density functional H(p) [2g, l27| . We have 
demonstrated here, however, that the integral in eq. (}4*Sj) can be exactly evaluated by the 
saddle point method provided that the thermodynamic limit is taken. The resulting real 
energy functional H(p) is given (for a homopolymer solution) by eq. ()35j) and involves the 
constrained (real) saddle point z0 o (r), which is the mean chemical potential field required 
to generate a particular density pattern p{r). 

We believe that our discovery will have a number of significant implications. First of all, 
it provides theoretical justification for ad hoc approaches that have been used to include 



24, 



281 ] . Secondly, by 



thermal fluctuations in density field-based simulations of polymers 
circumventing the sign problem, it avoids the difficulties that have been encountered in field- 
theoretic polymer simulations based on eq. (J4~§|) . Moreover, with a real field theory, a much 
wider variety of options exist for carrying out stochastic computer simulations. Finally, we 
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note that our proof that the <p integral in eq. ()48|) is dominated for large systems by the 
saddle point carries over to a much broader class of classical fluids, both simple and complex. 
For example, real density functionals can be constructed using the present methods for 
alternative polymer chain models, such as the worm-like chain used to describe semiflexible 
polymers. Real density-based field theories for a wide variety of polymers, copolymers, 
and alloys can be derived, including solutions of charged polyelectrolyte systems. In the 
case of simple fluids, we believe that the present formalism will provide new insights and 
numerical strategies for improving density functional methods, including the treatment of 
crystallization and melting. 

A clear disadvantage of the density field approach, at least for polymers, is that force or 
energy evaluations require the computation of the mean- field i(f> (r) from equations such as 
(JTUJl at constrained monomer density p(r). This must be done iteratively and each iteration 
requires an independent solution of the diffusion equation (JTTj). In constrast, evaluation of 
the energy Q{4>) in the potential field formulation of eq. ()49)1 requires only a single solution of 
the diffusion equation. Since the vast majority of the computational effort in a field-theoretic 
computer simulation of a polymeric fluid is spent in solving the diffusion equation, it is 
essential that a very efficient scheme be devised for computing i(p given p. We are cautiously 
optimistic, however, that the extra computational burden associated with simulating in the 
density variables will be more than offset by the avoidance of the sign problem. 
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